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(Dated: May 29, 2007) 



We present a measurement of the top quark mass in the dilepton channel based on approximately 
370 pb~ x of data collected by the D0 experiment during Run II of the Fermilab Tevatron collider. 
We employ two different methods to extract the top quark mass. We show that both methods yield 
consistent results using ensemble tests of events generated with the D0 Monte Carlo simulation. 
We combine the results from the two methods to obtain a top quark mass mt = 178.1 ± 8.2 GeV. 
The statistical uncertainty is 6.7 GeV and the systematic uncertainty is 4.8 GeV. 
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PACS numbers: 14.65.Ha 



The top quark mass is an important parameter in stan- 
dard model [H predictions. For example, loops involving 
top quarks provide the dominant radiative corrections to 
the value of the W boson mass. Precise measurements of 
the W boson and top quark masses provide a constraint 
on the Higgs boson mass [2j . 

At the Tevatron, top and antitop quarks are predom- 
inantly pair-produced. Top quarks decay to a W boson 
and a b quark. If the W bosons from the top and the anti- 
top quarks both decay leptonically (to ev or \xv) the final 
state consists of two charged leptons, missing transverse 
momentum (j/>t) from the undetected neutrinos, and two 
jets from the fragmentation of the b quarks. We call this 
the dilepton channel. It has a relatively small branching 
fraction (« 5%) but very low backgrounds. The mea- 
surement of the top quark mass in the dilepton channel 
is statistically limited. It provides an independent mea- 
surement of the top quark mass that can be compared 
with measurements in other tt decay channels, and a 
consistency check on the tt hypothesis in the dilepton 
channel. 

The D0 detector is a multipurpose collider detector [l| . 
The central tracker employs silicon microstrips close to 
the beam and concentric cylinders covered with scintil- 
lating fibers in a 2 T axial magnetic field. The liquid- 
argon/uranium calorimeter is divided into a central sec- 
tion covering pseudorapidity \ri\ < 1.1 and two endcap 
calorimeters extending coverage to \rj\ < 4.2 where 
r] = — m[tan(0/2)] and 6 is the polar angle with respect to 
the proton beam direction. The muon spectrometer con- 
sists of a layer of tracking detectors and scintillation trig- 
ger counters between the calorimeter and 1.8 T toroidal 
iron magnets, followed by two similar layers outside the 
toroids. 

We present here two measurements that were carried 
out independently by two groups of analyzers. Both 
groups chose to optimize their analyses in different ways, 
one using a relatively loose event selection, the other tak- 
ing advantage of the low background in top-antitop sam- 
ples selected using tagging of 6-quark jets. In the end, 
we combine the results from both analyses taking into 
account the correlations between the results. 

The event selection is based on the measurement of the 
cross section for it-production in the dilepton channelQ 
with a few modifications. The analyses use about 370 
pb _1 of data from pp collisions at -^5=1. 96 TeV collected 
with the D0 detector at the Fermilab Tevatron collider. 

We select events with two oppositely charged, isolated 
leptons (e or /z) with transverse momentum px > 15 GeV 
and at least two jets with px > 20 GeV. Electron can- 
didates are isolated clusters of energy in the electro- 
magnetic section of the calorimeter that agree in their 
profile with that expected from electromagnetic show- 
ers, based on Monte Carlo simulations, and that are 
matched with a charged particle track reconstructed in 
the central tracker. Electrons must be either in the 



central calorimeter (pseudorapidity \rj\ < 1.1) or in the 
forward calorimeter (1.5 < \rj\ < 2.5). Muons are re- 
constructed as tracks in the muon spectrometer with 
\rj\ < 2, matched to a charged particle track in the cen- 
tral tracker. They must be isolated from other activity in 
the calorimeter and in the tracker. Jets are reconstructed 
with the improved legacy cone algorithm[6J with cone size 
AB, = y/Arj 2 + A<p 2 = 0.5 and are restricted to \rj\ < 2.5. 
All jets were corrected using the standard D0 jet energy 
scale corrections 7f. 

We distinguish e/i, ee, and \±\i events. For e/i events we 
require Ht > 122 GeV, where Ht is the scalar sum of the 
larger of the two lepton pt values and the px values of 
the leading two jets. For ee events we require sphericity 
S > 0.15 and missing transverse momentum fa > 35- 
40 GeV, depending on the dielectron invariant mass 
m(ee), and we reject events with 80 < m(ee) < 100 GeV 
to reduce the background from Z — > ee decays. Spheric- 
ity is defined as 1.5 times the sum of the first two eigenval- 
ues of the normalized momentum tensor calculated using 
all electrons, muons and jets in the event. 

For /i/i events we require inconsistency with the Z — > 
/z/i hypothesis based on the x 2 °f a kinematic fit. In 
some Z — > /i/i events a muon momentum is significantly 
mismeasured. These events are not consistent kinemati- 
cally with Z decays and they are therefore not eliminated 
by the kinematic fit. The mismeasured muon momen- 
tum gives rise to pt imbalance in the muon direction. 
We therefore require fa > 35 GeV if the azimuthal an- 
gle between the leading muon and the direction of fa, 
Acf)($T,n) < 175°. We tighten the fa requirement to 85 
GeV if the leading muon and the fa are approximately 
collinear in the transverse direction. 

For our mass measurements we use the following sam- 
ples of events. The "6-tag" sample consists of events that 
have at least one jet that contains a secondary vertex 
tag with transverse decay length significance A xy > 7 0] . 
This sample has very low backgrounds. The "no-tag" 
sample consists of events that have no such secondary 
vertex tags. The 26 events in these two samples consist 
of 20 e/i events, 5 ee, and 1 /i/i event. 

The "tight" sample does not use the 6-tagging infor- 
mation. It contains all ee and /i/i events that are in either 
the 6-tag or the no-tag samples. For e/i events the tight 
sample requires the more restrictive cuts Ht > 140 GeV, 
fa > 25 GeV and tighter electron identification cuts 
to reduce backgrounds. To increase the acceptance for 
dilepton decays, we also analyze a sample of events that 
requires only one well-identified lepton (e or /i) with 
Pt > 15 GeV and an isolated track with pt > 15 GeV 
instead of the second identified lepton. The events must 
also have at least one jet with a secondary vertex tag, and 
fa > 15-35 GeV, depending on lepton flavor and the in- 
variant mass of the lepton+track system. We call this 
the "£+track" sample. Events with two well-identified 
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leptons are vetoed from this sample so that there is no 
overlap between the ^+track sample and the other dilep- 
ton samples. There are 6 e+track events and 3 /i+track 
events in this sample. The observed event yields for each 
of the data samples are listed in Table HI 

TABLE I: Expected and observed dilepton event yields for tt 
production with mt = 175 GeV and the backgrounds from 
WW and Z production based on Monte Carlo, and from 
misidentified leptons (mis-id) based on collider data. 



Sample 


ft 


WW 


Z 


Mis-id 


Total 


Data 


U no-tag 


7.2 


1.1 


2.6 


2.2 


13.2(1-) 


12 


U 6-tag 


9.9 


0.05 


0.12 


0.09 


10.1 ±0.9 


14 


11 tight 


15.8 


1.1 


2.4 


0.5 


19.8 ±0.6 


21 


^±track 


6.3 


0.01 


1.8 


0.4 


8.5 ±0.3 
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Monte Carlo samples are generated for nineteen values 
of the top quark mass between 120 and 230 GeV. The 
simulation uses ALPGEN [9( with CTEQ5L parton distri- 
bution functions as the event generator, PYTHIA [llj 
for fragmentation and decay, and GEANT [12| for the de- 
tector simulation. No parton-shower matching algorithm 
was used in the generation of these event samples. We 
simulate diboson production with ALPGEN and pythia 
and Z/-f* — > tt processes with pythia. The number of 
expected events are determined by applying the selection 
cuts to these Monte Carlo event samples. These samples 
are corrected for lepton, jet and 6-tagging efficiencies de- 
termined from collider data. 

The tagging efficiency for 6-jets is measured in a data 
sample enhanced in heavy flavor jets by requiring at least 
one jet with a muon in each event. Monte-Carlo based 
corrections are applied to correct for sample biases. The 
probability to tag a light-flavor jet is measured from col- 
lider data using events with a secondary vertex with neg- 
ative decay length, meaning that the tracks forming the 
secondary vertex meet in the hemisphere that is on the 
opposite side of the primary vertex from the jet. 

The energy of Monte Carlo jets is increased by 3.4% in 
addition to the nominal jet energy scale corrections. This 
factor was determined by fitting the top mass and the jet 
energy scale in lepton±jets events and brings the invari- 
ant mass distribution of the two jets from the W boson 
decay in lepton±jets Monte Carlo events in agreement 
with that observed in the data. 

Event yield normalizations for Z — > ee and Z — > /i/i 
are obtained from data. The number of events with 
misidentified leptons is dominated by jets misidentified 
as electrons. We construct a likelihood discriminant to 
distinguish electrons from misidentified jets based on the 
shape of the energy cluster in the calorimeter and the on 
the matched track. We determine the contamination by 
misidentified jets in our sample by fitting the distribu- 
tion of this likelihood discriminant before we cut on it. 
Expected yields for signal and background are given in 
TableD 

We use only the two jets with the highest pr in this 



analysis. We assign these two jets to the b and b quarks 
from the decay of the t and t quarks. If we assume a 
value mt for the top quark mass, we can determine the 
pairs of t and t momenta that are consistent with the ob- 
served lepton and jet momenta and j£/>. A solution refers 
to a pair of top-antitop quark momenta that is consistent 
with the observed event. For each assignment of observed 
momenta to the final state particles and for each hypoth- 
esized value of mt, there may be up to four solutions. We 
assign a weight function w(m t ) to each solution, as de- 
scribed below. Events for which no solution exists are 
rejected from our data and Monte Carlo event samples. 
The event yields in Table [T] include this additional selec- 
tion requirement. Two events from the collider data are 
rejected with this requirement. 

We consider each of the two possible assignments of 
the two jets to the b and b quarks. We account for detec- 
tor resolutions by repeating the weight calculation with 
input values for the lepton and jet momenta that are 
drawn from the detector resolution functions for objects 
with the observed momenta. We refer to this procedure 
as resolution sampling. For each event we obtain a weight 

W(m t ) = 1/N x J2f=i Yh=\ w ( m t)i] by summing over 
all n solutions and averaging over N resolution samples. 
This weight characterizes the likelihood that the event is 
produced in the decay of a it pair as a function of m t . 

The techniques we use are similar to those used by the 
D0 Collaboration to measure the top quark mass in the 
dilepton channel using Run I data [13| . The data are 
analyzed using two different methods that differ in the 
event samples that they are based on, in the calculation 
of the event weight, and in the algorithm that compares 
the weights for the observed events to Monte Carlo pre- 
dictions to extract the top quark mass. 

The matrix-element weighting technique (A4WT) fol- 
lows the ideas proposed by Dalitz and Goldstein [14| and 
Kondo [l5[ . The solution weight is 

w(m t ) = f(x)f(x)p(E* t \m t )p{Ej\m t ), 

where f{x) is the parton distribution function of the pro- 
ton and x (x) is the momentum fraction carried by the 
initial (anti)quark. The quantity p(E^\m t ) is the prob- 
ability that the lepton has energy E% in the top quark 
rest frame for the hypothesized top quark mass mt- 

For each event we use the value of the hypothesized top 
quark mass m pea k at which W(m t ) reaches its maximum 
as the estimator for the mass of the top quark. We gener- 
ate probability density functions of m pca k for a range of 
top quark masses using Monte Carlo simulations. We call 
these distributions templates. To compute the contribu- 
tion of backgrounds to the templates, we use Z — > tt and 
WW Monte Carlo events. Backgrounds arising from de- 
tector signals that are misidentified as electrons or muons 
are estimated from collider data samples. 

We compare the distribution of rop ea k for the observed 
events to these templates using a binned maximum like- 
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lihood fit. The likelihood is calculated as 



L(m t ) = f[ 



i=l 



n s s % {m t ) + n b bi 



where rij is the number of data events observed in bin i, 
spirit) is the normalized signal template contents for bin 
i at top quark mass m t , bi is the normalized background 
template contents for bin i. The product runs over all 
rtbin bins. The background template consists of events 
from all background sources added in the expected rel- 
ative proportions. The signal-to-background fraction is 
fixed to n s /n b with the numbers of signal and background 
events (n a , rib) taken from Table |TJ 

To calibrate the performance of our method, we gener- 
ate a large number of simulated experiments for several 
input top quark mass values. We refer to each of these 
experiments as an ensemble. Each ensemble consists of as 
many events of each type as we have in our collider data 
sample. A given event is taken from the signal and back- 
ground samples with probabilities that correspond to the 
fraction of events expected from each sample. We use a 
quadratic function of mt to fit the — In L points to thir- 
teen mass points centered on the point with the smallest 
value of — InL. The distribution of measured top quark 
mass values from the ensemble fits gives an estimate of 
the parent distribution of our measurement. The ensem- 
ble test results indicate that the measured mass tracks 
the input mass with an offset of 1.9 ± 0.8 GeV, which we 
correct for in the final result. 

In general, the tails of the likelihood distribution for an 
ensemble are not well approximated by a Gaussian. Thus 
it is necessary to restrict the range of mass points that 
is included in the fit to points near the observed mini- 
mum in — InL. For small data samples, however, there 
is a substantial statistical uncertainty in the computed 
likelihood values which can be reduced by increasing the 
number of mass points used in the fit. Thus the range of 
mass points that are included in the likelihood fit must 
be optimized for the observed data sample size to obtain 
the best possible agreement between measured top quark 
mass and input top quark mass. This was done for both 
analyses based on Monte Carlo ensembles that contain 
exactly as many events as we observe in the data. 

The MWT analysis uses the no-tag and &-tag sam- 
ples of events. Separating out the very-low-background 
6-tagged events improves the precision of the result. The 
analysis is performed with separate templates for ee, e/x, 
and nn events and separate signal-to-background frac- 
tions for events without a 6-tag and > 1 6-tags. The 
maximum of the joint likelihood for all events, shown in 
Fig. [TJ corresponds to m t = 176.2 ± 9.2(stat) GeV after 
the offset correction. Figure [H shows the distribution of 
m pea k from collider data compared to the sum of Monte 
Carlo templates with mt — 180 GeV. 

The neutrino weighting technique (z/WT) ignores the 
measured $t in reconstructing the event. Instead we as- 
sume a representative range of values for the pseudora- 
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FIG. 1: Joint likelihoods from the MWT analysis (closed 
circles) and the i/WT analysis (open circles). The minima 
of the likelihood curves do not include the correction for the 
offset in the response. 
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FIG. 2: Distribution of m pca ^ from the .MWT analysis (cir- 
cles) compared to the sum of Monte Carlo templates for the 
no-tag and 6-tag channels and all lepton flavors for m t = 
180 GeV (open histogram). The shaded histogram indicates 
the background contribution. 
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pidities of the two neutrinos and the solution weight 



w(m t ) 



N, 



E 



exp 



■n i=i 



2al 



exp 



-j>v? 



2o* 



characterizes the consistency of the resulting solutions 
with the observed $t- The sum is over the N v steps of 
neutrino rapidity values, $ x , andj^ are the x and y com- 
ponents of the sum of the neutrino momenta computed 
for step i, and a x and <r y are the measurement resolu- 
tions for ]/> x andj^j,. We then normalize the event weight 
W(m t ) over the range 80 < m t < 330 GeV and integrate 
it over ten bins in m t . Every event is thus characterized 
by a 9-component vector W = (Wi, Wg) (the 10 th bin 
is fixed by the first nine and the normalization condition). 
We compare the vectors from the collider data events to 
sets of N Monte Carlo events generated with different 
values of mt by computing the signal probability 



N 9 



f s (W\m t ) = -Y,T] 



e^[-{W l -W l f c f/2h 2 ] 
* h l\ Jo exp[-(T4^' - W!? c )*/2h*]dW> ' 

where W^ 1C is the vector of weights from Monte Carlo 
event j. The value of the resolution parameter h is op- 
timized using ensemble tests based on simulated events 
to give the best agreement between input mass and mea- 
sured mass. We compute a similar probability fb(W) for 
backgrounds and combine them in the likelihood 

L(m t ,nb, n) = G(rib - rib, a)P{n s + rib, n) 
n s f s (Wi\m t ) +n b fb(Wi) 



n 



rib 



which we optimize with respect to mt, the number of 
signal events n s , and the number of background events 
rib- G is a gaussian constraint on the difference between 
rib and the expected number of background events n b , 
and P is a Poisson constraint on n s + rib to the number 
of events n observed in data. 

The ^WT analysis uses the tight sample and the 
£+track sample. The analysis is performed with separate 
templates for ee, efx, and [ill events in the tight sample 
and the two lepton flavors in the £+track sample. We fit 
the — In L points for values of within 20 GeV of the 
point with the smallest value of — In L with a quadratic 
function of m t . The performance of the vWT algorithm is 
checked using ensemble tests as described for the MWT 
algorithm. The average measured values of mt track the 
input values with an offset of 1.7±0.2 GeV. For the z/WT 
analysis, the maximum of the joint likelihood of all events 
(Fig. [JJ) corresponds to m t = 179.5 ± 7.4(stat) GeV after 
the offset correction. 

We also use ensemble tests to study the size of sys- 
tematic uncertainties (see Table [TTJ) . By far the largest 
systematic uncertainty originates from the uncertainty in 
the calibration of the jet energy scale of 4.1%. We deter- 
mine the effect of the uncertainty on the measurement 



by generating ensemble tests with the jet energy scale 
increased and decreased by one standard deviation. 

We estimate the sensitivity of the result to uncertain- 
ties in the parton distribution functions by analyzing en- 
sembles generated with a range of available parton dis- 
tribution function sets. The next to largest uncertainty 
originates from the modeling of gluon radiation in the 
Monte Carlo. Gluon radiation can give rise to additional 
jets in the event. In the data about one third of the events 
have more than two jets. The two analyses used different 
procedures to estimate this effect. For the ^WT analysis, 
events with three reconstructed jets from tt + 1 parton 
events generated with ALPGEN were analyzed in ensemble 
tests with templates derived from tt events with only two 
jets and the difference in reconstructed top quark mass 
was applied as an uncertainty to the fraction of events 
with more than two jets. In the A4WT analysis the frac- 
tion of events with only two jets was varied in ensem- 
ble tests within the range that is consistent with the jet 
multiplicity spectrum observed in the data and analyzed 
with the nominal templates. The observed variation in 
the result was applied as systematic error. 

We estimate the effect of uncertainties in the shape 
of the background distributions to determine the back- 
ground uncertainty. For the AiWT analysis we also per- 
form tests with ensembles in which we varied the back- 
ground fraction, which was fixed in the mass fit, by its 
uncertainty. For the -^-l-track sample, the heavy flavor 
content in the background is a significant source of un- 
certainty. This only contributed to the vWT analysis. 
The finite size of the Monte Carlo samples limits the sta- 
tistical precision with which we can extract the top quark 
mass. This is accounted for in the Monte Carlo statistics 
uncertainty. Finally, we generated ensembles with var- 
ied jet and muon momentum resolutions to estimate the 
effect of their uncertainties. The resulting uncertainties 
for the i/WT analysis are quoted in Table [TTJ The effect 
on the result of the VWWT analysis was negligible. 

We follow the method for combining correlated mea- 
surements from Ref. [HI ] in combining the results from 
the AiWT and i/WT analyses. We determine the sta- 
tistical correlation between the two measurements using 
ensemble tests. The correlation factor between the two 
analyses is 0.35. The systematic uncertainties from each 
source in Table [TT] are taken to be completely correlated 
between the two analyses. The results of the combination 
are also listed in Table UTl 

In conclusion, we measure the top quark mass in the 
dilepton channel. We obtain m t = 178.1 ± 6.7(stat) ± 
4.8(syst) GeV as our best estimate of the top quark mass. 
This is in good a gree ment with the world average m t — 
172.5 ± 2.3 GeV 0, based on Run I and Run II data 
collected by the CDF and D0 Collaborations. 

We thank the staffs at Fermilab and collaborating in- 
stitutions, and acknowledge support from the DOE and 
NSF (USA); CEA and CNRS/IN2P3 (France); FASI, 
Rosatom and RFBR (Russia); CAPES, CNPq, FAPERJ, 
FAPESP and FUNDUNESP (Brazil); DAE and DST 



8 



TABLE II: Summary of dilepton mass measurements. 





MWT 


uWT 


Combined 


Top quark mass 


176.2 


179.5 


178.1 


GeV 


Statistical uncertainty 


9.2 


7.4 


6.7 


GeV 


Systematic uncertainty 


3.9 


5.6 


4.8 


GeV 


Jet energy scale 


3.6 


4.8 


4.3 


GeV 


Parton distribution functions 


0.9 


0.7 


0.8 


GeV 


Gluon radiation 


0.8 


2.0 


1.5 


GeV 


Background 


0.2 


1.4 


0.9 


GeV 


Heavy flavor content 




0.6 


0.3 


GeV 


Monte Carlo statistics 


0.8 


1.0 


0.9 


GeV 


Jet resolution 




0.6 


0.3 


GeV 


Muon resolution 




0.4 


0.2 


GeV 


Total uncertainty 


10.0 


9.3 


8.2 


GeV 
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